Extraordinary electron transmission through a periodic array of quantum dots 
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We study electron transmission through a periodic array of quantum dots (QD) sandwiched 
between doped semiconductor leads. When the Fermi wavelength of tunneling electron exceeds 
the array lattice constant, the off-resonant per QD conductance is enhanced by several orders of 
magnitude relative to the single-QD conductance. The physical mechanism of the enhancement is 
derealization of a small fraction of system eigenstates caused by coherent coupling of QDs via the 
electron continuum in the leads. 



Interference effects in electron transmission through 
localized states such as, e.g., semiconductor quantum 
dots (QD) have been among the highlights in coherent 
transport studies [l|, 0] ■ The phase acquired by an elec- 
tron tunneling along multiple pathways provided by QDs 
placed between doped semiconductor leads can cause 
striking features in tunneling current Q. The simplest 
realization of coherent transport is provided by two QDs 
independently coupled to the leads at weak (or absent) 
interdot tunneling [J]. The conductance of a double-QD 
system in a magnetic field exhibits Aharonov-Bohm os- 
cillations [5|-[7| as a function of magnetic flux penetrating 
the area bound by tunneling paths @, H, 0] . At zero field, 
the coherence is controlled by QDs' coupling via the elec- 
tron continuum in the leads [4.ll0}. If QDs separation, a, 
is comparable to electron's Fermi wavelength, Xp, then 
electron transmission is mediated by the system eigen- 
states rather than individual QDs, leading to conduc- 
tance peak narrowing or Fano-like lineshapes [i, 11- 13 1 . 
A close optical analogy for coherent transport through 
QDs is the Dicke superradiance of two excited atoms [14j ; 
QDs' coupling via the electron continuum is similar to ra- 
diative coupling between the atoms 0, E3] • 

In this Letter, we study electron transport through 
a periodic array of localized states, e.g., a QD lattice 
with period a separated from the electron gas (EG) in 
the leads by tunneling barriers. We consider a linear 
chain of QDs situated between two-dimensional (2D) 
leads (1DQD-2DEG) and a square lattice of QDs sand- 
wiched between 3D leads (2DQD-3DEG), as illustrated 
in Fig. [TJ We demonstrate that in the subwavelength 
regime, i.e., a/Xp < 1, the off-resonant per QD conduc- 
tance is strongly enhanced relative to single-QD one. We 
trace the origin of this effect to delocalization of a small 
fraction of system's eigenstates which mediate electron 
transport across tunneling barriers. 

The novel coherent transport effect predicted here 
bears similarity to extraordinary optical transmission 
(EOT) of light through a metal film perforated with a 
periodic array of nanoholes [15j . Suggested EOT mech- 
anisms involve surface plasmon-polaritons TB-T3|, dy- 
namical diffraction [3], light's coupling to waveguide res- 
onances [III], and composite diffracted evanescent wave 
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FIG. 1. (Color online) Schematics of (a) ID periodic chain of 
QDs between 2D leads (1DQD-2DEG), (b) 2D lattice of QDs 
between 3D leads (2DQD-3DEG), and (c) lowest EG bands 
(g — 0, ±27r/a) in k z — plane in units of E a — 7r 2 /2ma 2 ; 
vHS is located at the crossings of bands and Fermi level. 



2fJ. However, the physic al origin of EOT remains under 



active discussion [21| - |23| . A similar effect was recently 
reported for a transmission of sound waves 24- 2(|. Here 
we describe a related phenomenon in electron transport 
through a QD array. Despite differences with EOT, the 
generic effect of enhanced transmission persists. In this 
case, however, we can pinpoint the mechanism of extraor- 
dinary electron transmission (EET). 

Indeed, consider the electronic states in the system 
comprised of QD lattice (ID or 2D) tunnel-coupled to 
EGs in the leads (2D or 3D) without direct tunneling 
between QDs (see Fig. [I]). The periodic potential of QD 
lattice splits the EG momentum space (k,k z ) into bands 
in transverse direction (p + g,k z ), where p lies in the 
first Brillouin zone (BZ), g are reciprocal lattice vectors, 
and k z is the momentum along tunneling direction. At 
Fermi level, the electronic states in band g with energy 
Ep = (p + g) 2 /2m + k z /2m are either propagating for 
\p + g\ < kp (real k z ), or evanescent for \p + g\ > kp 
(imaginary k z ), where kp — \J2mEp = 2tt/Xf and m are 
the electron Fermi momentum and mass, respectively, 
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(we set h = 1). At the same time, the coupling be- 
tween QDs via EG gives rise to localized states' band 
(LSB) whose eigenstates are characterized by transverse 
quantum number p and are confined along z direction. 
The LSB eigenstates' dispersion E p comes from coher- 
ent coupling between QDs' via evanescent states of EG, 
while eigenstates' decay into EG propagating states de- 
termines their lifetime. Importantly the conservation of 
p reduces the phase space for tunneling to that of a dis- 
crete set of IDEGs (enumerated by g) with dispersion 
k 2 /2m. A transition between an LSB eigenstate p and 
EG is then determined by the 1DEG partial density of 
states (DOS) which has the van Hove singularity (vHS) 
at \k z \ = \fkp — (p + g) 2 = (see Fig. [T|). An electron 
in LSB eigenstate with p near vHS can escape to the 
leads and hence is delocalized across tunneling barriers. 

Specifically a single-QD resonant tunneling conduc- 
tance through level Eq is given by standard expression 
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where Tl,r are decay rates to left or right leads. In the 
symmetric case Tl = L/? = T, Go = e 2 /w at resonance 
Ef = Eq, while the off-resonant Go is reduced by factor 
{T/SE) 2 < 1, where 6E — Ep — E is the Fermi level 
detuning. We show that the low-energy per QD array 
conductance G is enhanced relative to Gq by the factor 
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7 = (27r /akp)T being fop-independent, which can reach 
several orders of magnitude. Equation © holds in the 
2DQD-3DEG case; in the 1DQD-2DEG case the prefac- 
tor is replaced by A/nakp and 7 = 2r. In both cases, 
the prefactor reflects QDs' number in the volume with 
linear size Xf and results from the coherent coupling be- 
tween QDs 27|, while the logarithmic factor describes 
delocalized states' contribution. For asymmetric leads, 
enhancement is reduced due to different vHS locations in 
p-space for left and right leads. 

Theory. — We start with tunneling Hamiltonian for a 
square lattice of N QDs with in-plane coordinates tv, sep- 
arated by potential barriers from 3DEGs in the leads 



H = Eoc)cj + £?4«c„ Q + Y (yfrcfaa + H.c 
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where 



(3) 

and Eq are, respectively, creation and an- 
nihilation operators and energies for localized states, 
ci a , c vai and £" are those for EG states in the leads 
(a = L, R), and V„" is transition matrix element between 
localized and EG states. We assume the barrier thickness 
smaller than the lattice period and neglect the direct tun- 
neling between QDs. The zero-temperature conductance 



normalized per QD is given by [lfj| 

G = ff * f L 

irh N \ Ef-Eq-Y, Ep-Eq-T,^ 

where = Ef- + E^ is the self-energy matrix, 
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whose principal and singular parts determine the energy 
matrix Ay and the decay matrix rg., respectively, and 
the trace is taken over lattice sites. The matrix element 
can be presented as V°" v = (AL z )~ 1 / 2 e lk " r H a , where fe 
and rj are the in-plane (transverse) momentum and co- 
ordinate, respectively, t a is the tunneling amplitude be- 
tween QD and the lead a (we assume that the barrier is 
sufficiently high and neglect its dependence on energy), 
and A = Na 2 and L z are normalization area and length, 
respectively. The self-energy is expressed via the electron 
Green function as = £„G a (rj, d a ; rj, d a ), d a being 
the barrier thickness. 

The coupling between QD lattice and EG gives rise 
to LSB whose eigenstates are characterized by in-plane 
quasimomentum p lying in the first BZ. At the same time, 
QDs' periodic potential splits 3DEG momentum space 
(fe, k z ) into bands (p+g, k z ), where g = (2wm/a, 2nn/a), 
to and 77 being integers. The LSB spectrum is obtained 
via Fourier transform of Eq. ([5]) , 
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[Ep - (k 2 + k 2 z )/2m a ] 1 are, re- 
spectively, the LSB self-energy and the Fourier of elec- 
tron Green function, m a being the electron mass. After 
fez-integration, the LSB self-energy takes the form 
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where j a = 2m a t 2 a /a, and kp = y/2m a Ep. Note that 
for a given p, each band contributes either to Ap or to 
Tp depending on whether |gr +p\ > k F or \g + p\ < kp. 
Thus, QDs' coupling via evanescent states of EG gives 
rise to the LSB dispersion E p = E + A p + A p (note that 
A^ < 0), while transitions between LSB and propagating 
states of EG determine the decay rate . The LSB self- 
energy Eq. ([7]) diverges at |fl+p| = k F , i.e., for EG states 
propagating along the interface, due to vHS in the 1DEG 
DOS at k z = 0. The latter is caused by conservation of 
p which reduces the phase space for tunneling to that of 
1DEG with dispersion k 2 /2m. An LSB eigenstate with 
p near vHS has a high decay rate to the EG and is, in 
effect, delocalized across the barriers. As we show below, 
these states mediate EET through a QD lattice. 
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Using Fourier transform ([6]), the conductance (|4|) can 
be recast in terms of system's eigenstates as 
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where p-integral is taken over 2D BZ. Thus, G is a sum 
of partial conductances for each p: G — JdpG p , 
where G p has the form ([T]) with Eq and T a replaced by 
E p and Tp, respectively. In 1DQD-2DEG case, p and 
g in Eq. (0 are ID vectors (and j a = 2m Q t^), while 
integration measure in Eq. © extends over ID BZ. 

For a given Ep, only EG bands with propagating 
states, i.e., those with \g + p\ < k F , contribute to Tp 
and hence provide transmission channels. If Ep is larger 
than the geometric energy scale E a — ir 2 /2ma 2 , i.e., the 
coherent coupling between QDs is weak (akp 3> 1), then 
a large number of bands contribute to G, and replac- 
ing the sum over g in Eq. (0 with integral one recovers 
the single-QD decay rate Tp = r Q and conductance ([1]) 
for uncoupled QDs. However, for Ep < E a , i.e., strong 
coherent coupling (akp < 1), only lowest EG bands con- 
tribute to decay rate T^, and G is determined by its 
behavior near vHS. 

In Fig. [2] we show the conductance and enhancement 
factor for symmetric 1DQD-2DEG case (Tp = T p , j a — 
7, k F = kp). Numerical calculations were carried using 
Eqs. ([7]) and (|S} with QD lattice period chosen to set 
Eq = 9E a , corresponding io akp = 3-7T at resonance, and 
g with m,n up to 50 included. For akp /it < 1, only 
g = band contributes to T p ; T p ^ only below vHS 
(p < kp), while E p is affected by vHS for p > kp [see 
Fig. EJc)]. Note that higher bands' contribution leads 
to energy renormalization by a constant that causes ap- 
parent deviation of E p /Eq from unity below vHS (for 
the given maximal g); in the following, we absorb this 
constant into Eq. Note also that E p /Eq is suppressed 
relative to F p /r by small factor T/Eq -C 1, and so the 
LSB dispersion E p has no significant effect on G. In 
fact, the LSB DOS (not shown here) is nearly indistin- 
guishable from that of a single QD, indicating that vHS 
affects only a small fraction of LSB. However, it is this 
small fraction of LSB eigenstates p with large decay rates, 
-Tp ^ 1*^1 > that causes dramatic conductance enhance- 
ment on the low energy side [see Fig. [Ha)]. Indeed, the 
partial conductance of these off-resonant states is close 
to the maximal (ballistic) value [G p ~ e 2 /7r in Eq. ([8]l] 
indicating that these states are, in effect, delocalized. For 
akp /it < 1, G can be evaluated as (in 2DQD-3DEG case) 
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where chief contribution to the integral comes from the 
region (p — kp)/kp ~ (7 / 'SEakp) 2 <C 1 constituting 




FIG. 2. (Color online) 1DQD-2DEG normalized conductance 
(a) and enhancement factor (b) vs. EG Fermi momentum are 
shown for symmetric case at Eo/E a ~ 9.0 for two values of 
7/-E0; (c) vHS in LSB energy spectrum at a/cp/V = 0.2. 



~ (-y/SE) fraction of 2D BZ, while the log factor comes 
from the slow (square root) divergence of vHS and re- 
flects delocalized states' distribution in energy. In 1DQD- 
2DEG case, G differs from the result © by factor 2/ akp. 
The sharp features in fc^-dependence of G occur when 
the Fermi level crosses EG band edges. The enhance- 
ment factor G/Gq is now determined by taking into 
account the relations 7 = 2r for 1DQD-2DEG and 
7 = (27r/afc F )r for 2DQD-3DEG case, yielding G/G Q = 
(4/irakp) In (akp \ SE\/y) in the former case and Eq. @ 
in the latter case. In the region j/\5E\ <C kpa <C 1, 
both logarithm and its prefactor contribute to overall en- 
hancement that reaches 2 orders of magnitude in 1DQD- 
2DEG case [see Fig.[^b)]; however, the (akp)~ 2 prefac- 
tor dependence in Eq. © leads to greater enhancement 
in 2DQD-3DEG case, as shown in Figs.[3ja) and[3jb). 

Breaking the left-right symmetry significantly alters 
the conductance; however, a distinction must be made 
in the way it is broken [see Fig. [3tc) and EJd)]. For 
example, a smaller Tp due to higher tunneling barrier 
between QDs and the right lead causes overall reduction 
of both array and single-QD conductances; however the 
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FIG. 3. (Color online) 2DQD-3DEG enhancement factor (a) 
and normalized conductance (b) vs. EG Fermi momentum 
are shown for symmetric case at Eo/E a — 9.0 for two val- 
ues of 7/-E0; (c) and (d) - symmetric case is compared to 
asymmetric case for jl/Eq = 10~ 2 . 



enhancement factor G/Gq stays essentially unchanged. 
In contrast, an EG asymmetry, e.g., due to different ef- 
fective masses (mi 7^ wir), drastically lowers the en- 
hancement. The reason is that for kp 7^ kp, the LSB 
eigenstates have different decay rates to left and right 
EGs due to different respective vHS locations in the p- 
space. Here the enhancement factor is easily evaluated as 
G/G « 2ir/(ak F ) 2 v / m R /m L - 1 for (m R /m L - 1) ~ 1, 
i.e., it is suppressed by log factor relative to Eq. 

Finally, consider the effect of a finite temperature, T, 
and of a weak disorder due to QDs' size and position 
distributions. The disorder can be characterized by in- 
homogeneous width of QDs energies, 70 <C 5E, and devi- 
ation of lattice constant from its average value, Sa <C a. 
The former lowers the resonant transmission but does 
not affect the off-resonant one dominated by states with 
T p SE. Similarly, finite T broadens the resonance, but 
the states with r p 3> T remain unaffected. However, 
even a weak nonperiodicity breaks down exact conserva- 



tion of p and hence smears out vHS thus reducing the 
fraction of delocalized states. It is easy to show that for 
<W a £ il/ SE f, E q- © still holds , but with the loga- 
rithm argument replaced by y a/ 5a. 

In conclusion, there is an apparent similarity between 
EET through a QD lattice and EOT picture based on 
the couplin g of light to waveguide modes in a perforated 
metal film [19]. Note, however, that despite strong en- 
hancement (up to I0 4 in Fig. [3]), the off-resonant conduc- 
tivity is still much lower than its resonance value, in con- 
trast to EOT. The reason can be traced to exponentially 
small barriers' transparencies which render the fraction 
of delocalized states to be exponentially small as well. 
On the other hand, EOT represents a series of well de- 
fined peaks while EET through a QD lattice reveals itself 
in elevated plateau on the conductance low-energy side. 
These and other issues will be addressed by us elsewhere. 
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